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We present a search for the decay — > t~^v using 383 x 10^ BB pairs collected at the T{AS) 
resonance with the BABL4R detector at the SLAC PEP-II B-Factory. A sample of events with one 
reconstructed semileptonic B decay {B~ D°£~PeX) is selected, and in the recoil a search for 

B+ ^ r+ V is performed. The r is identified in the following channels: t"*" — » e^vv, r+ [i^vv, 
r"*" — > T^^V and n'^n'^V. We measure a branching fraction of B{B^ — > t'^v) = (0.9±0.6(stat.)± 



4 



0.1(syst.)) X 10 In the absence of a significant signal, we calculate an upper limit at the 90% 
confidence level of B{B'^ t^v) < 1.7 x 10^''. We calculate the product of the B meson decay 
constant /s and \Vub\ to be /s • \Vub\ = (7.2l^;°(stat.) ± 0.2(syst.)) x 10"* GeV. 

PACS numbers: 13.20.-v, 13.25.Hw 



INTRODUCTION 



II. THE BABAR DETECTOR AND DATASET 



In the Standard Model (SM), the purely leptonic decay 
^ r+J^ [1] proceeds via quark annihilation into a 
boson. The branching fraction is given by: 



Stt 



rrir 



(1) 

where Vub is an element of the Cabibbo-Kobayashi- 
Maskawa quark- mixing matrix [2, 3], /s is the B me- 
son decay constant, Gi? is the Fermi constant, tb+ is the 
lifetime, and and rrir are the 5+ and r masses. 
Physics beyond the SM, such as two-Higgs doublet mod- 
els, could enhance or suppress B{B~^ t^u) through 
the introduction of a charged Higgs boson [4-6]. Using 
theoretical calculations of fs from lattice QCD and ex- 
perimental measurements of \Vub\ from scmilcptonic B 
decays, this purely leptonic B decay can be used to con- 
strain the parameters of theories beyond the SM. Or, as- 
suming that SM processes dominate and using the value 
of \Vub\ determined from semileptonic B decays, purely 
leptonic decays provide an experimental method of mea- 
suring /s with reduced theoretical error. 

The branching fractions for B^ ^i^v and i?+ e^v 
are suppressed by factors of --^ 5 x 10~^ and 10"'^ 
with respect to B+ — > t^u. However, a search for 
5+ r+z^ is experimentally challenging due to the 
large missing momentum from multiple neutrinos, which 
makes the signature less distinctive than in the other lep- 
tonic modes. The SM estimate of the branching fraction 
for B+ T+v, using l^fcl = (4.31 ± 0.30) x 10-^ [7] 
and /b = 0.216 ± 0.022 GeV [8] in Eq. 1 is (1.6 ± 
0.4) X 10"''. In a previously published analysis using 
a sample of 223 x 10^ r(45) decays, the BABm collab- 
oration set an upper limit of B{B'^ t^v) < 2.6 x 
10-* at the 90% confidence level (CL) [9]. The Belle 
Collaboration has reported evidence from a search for 
this channel where the branching fraction was measured 
to be B{B+ T+v) = (1.79tg:^^(stat.);°;^^(syst.)) x 
10-4 [10]. 



* Deceased 

^Also witfi Universita di Perugia, Dipartimento di Fisica, Perugia, 
Italy 

tAlso witfi Universita dclla Basilicata, Potenza, Italy 

^Also with Univcrsitat dc Barcelona, Facultat de Fisica, Departa- 

incnt ECM, E-08028 Barcelona, Spain 

^Also with IPPP, Physics Department, Durham University, 
Durham DHl 3LE, United Kingdom 



The data used in this analysis were collected with the 
BAfilK, (k^tcctor [11] at the PEP-II storage ring. The sam- 
ple corresponds to an integrated luminosity of 346 fb^^ at 
the T(45) resonance (on-resonance) and 36.3 fb"^ taken 
at 40MeV below the BB production threshold (off- 
resonance) which is used to study background from 
e+e" ff if = u,d,s,c,T) continuum events. The 
on-resonance sample contains (383 ± 4) x 10^ T{4S) de- 
cays. The detector components used in this analysis are 
the tracking system composed of a five-layer silicon ver- 
tex detector and a 40-layer drift chamber (DCH), the 
Cherenkov detector for charged w-K discrimination, a 
Csl calorimeter (EMC) for photon and electron identifi- 
cation, and an 18-layer flux return (IFR) located outside 
of the 1.5 T solenoidal coil and instrumented with resis- 
tive plate chambers for muon and neutral hadron identi- 
fication. For the most recent 133 fb~^ of data, a portion 
of the resistive plate chambers has been replaced with 
limited streamer tubes. We analyze the data from sev- 
eral data-taking periods separately to account for varying 
accelerator and detector conditions. 

A GEANT4-based [12] Monte Carlo (MC) simulation 
is used to model signal efficiencies and physics back- 
gromids. The t Icpton decay is modeled using EvtGen 
[13]. Beam-related background and detector noise from 
data are overlaid on the simulated events. Simulation 
samples equivalent to approximately three times the ac- 
cumulated data are used to model BB events, and sam- 
ples equivalent to approximately 1.5 times the accumu- 
lated data are used to model continuum events. We de- 
termine selection efficiencies for signal events using a MC 
simulation where one B~^ meson decays to t^v, while the 
other is allowed to decay into any final state. 



III. ANALYSIS METHOD 

Due to the presence of miiltiplc neutrinos, the B"^ 
T'^v decay mode lacks the kinematic constraints which 
are usually exploited in B decay searches in order to re- 
ject both contimmm and BB backgrounds. The strategy 
adopted for this analysis is to reconstruct exclusively the 
decay of one of the B mesons in the event, referred to as 
the "tag" B. The remaining particle(s) in the event (the 
"recoil" ) are asstuned to come from the other B and are 
compared with the signature expected for — > t'^v. In 
order to avoid experimenter bias, the signal region in data 
is blinded until the final yield extraction is performed. 

The tag B is reconstructed in the set of scmilcptonic 
B decay modes B~ ^ D^£~i>^X, where £ denotes either 
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electron or muon, and X can be cither nothing or a tran- 
sition particle (tt^ or photon) from a higher mass charm 
state decay which we do not attempt to explicitly include 
in the tag B. However, we explicitly veto events where 
the best tag candidate is consistent with neutral B decay, 
where the X system is a single charged pion that can be 
combined with the -D° to form a £>*+ candidate. 

The T'^v signal is searched for in both leptonic 

and hadronic r decay modes constituting approximately 
71% of the total r decay width: t+ e^vV, r"*" — » n'^vV, 
r+ — > TT+F, and r+ — > Tr+TT^F. We do not consider the 
7r+7r~7r+i7 mode since we found it to be dominated 



xlO' 



xlO' 



by background events. 



A. Tag B Reconstruction 

D'^i candidates are reconstructed by combining a 
with an identified electron or muon with momentum 
above 0.8GeV/c in the e+e~ center-of-mass (CM) frame 
(Fig. 1). The flight direction of the is required to 
intersect with the lepton track. Assuming that the mass- 
less nciitrino is the only missing particle, we calculate the 
cosine of the angle between the D^i candidate and the 
B meson, 



cos^ 



B-DOe 



m. 



m 



2\p 



b\\Pd» 



(2) 
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FIG. 1: D" mass for tag B candidates containing an (a) elec- 
tron or (b) muon and the CM momentum of the tag B lepton 
for tag B candidates containing an (c) electron or (d) muon. 
On-resonancc data (filled circles) arc overlaid on the BB MC 
(solid histogram) and non-resonant background MC (gray his- 
togram), which have been normalized to the integrated data 
luminosity. Off-resonance data (open diamonds) are over- 
laid for comparison, and normalized to the on-resonance in- 
tegrated luminosity. 



where {Ejjoi, pnoe) and {Eb, Pb) are the four-momenta 
of the D^(. and B in the CM frame, and m^jo^ and 
TTiB are the masses of the D^t candidate and B'^ me- 
son (the nominal mass [7] is used), respectively. Eb and 
the magnitude of pb are calculated from the beam en- 
ergy: Eb — -Ecm/2, where i?cM is the CM energy of the 
beams, and \pb\ 



Correctly reconstructed 



candidates populate cosOb-doi in the range of [—1,1], 
whereas combinatorial backgrounds can take unphysical 
values outside this range. We retain events in the interval 
—2.0 < cosOB-D°e < l-l; where the upper bound takes 
into account the detector resolution and the less restric- 
tive lower bound accepts those events where the X is a 
soft transition particle from a higher mass charm state. 
We reconstruct the candidates in four decay modes: 



K- 



K- 



K-TT+TT^. and K%+-k- 



only con- 



sidering K^: candidates decaying to charged pions. The 
charged tracks are required to meet particle identifica- 
tion criteria consistent with the particle hypothesis and 
are required to converge at a common vertex. The 7r° 
candidates are required to have invariant masses between 
0.115 and 0.150 GeV/c^ and the photon daughter candi- 
dates of the 7r° must have a minimum laboratory en- 
ergy of 30 MeV and have shower shapes consistent with 
electromagnetic showers. The mass of the reconstructed 
candidates (Fig. 1) in the K~it^ , K~n^Tr~n^ , and 
ii'gTr+Tr" modes is required to be within 20 MeV/c^ of the 
nominal mass [7], while in the K~i:~^'n:^ decay mode the 
mass is required to be within 35 MeV/ of the nominal 



mass. These constraints are determined by fitting a sin- 
gle Gaussian function and a first-order polynomial to the 
mass distribution in signal MC and correspond to the Scr 
positions on the Gaussian. Furthermore, the sum of the 
charges of all the particles in the event must be equal 
to zero. If more than one suitable D^i candidate can 
be reconstructed, the best candidate is taken to be the 
one with the largest probability of converging at a single 
vertex. 

The tag reconstruction efficiency, extracted from sig- 
nal MC and averaged over all data taking periods, is 
(6.64 ± 0.03) X 10~^, where the error is due to the statis- 
tics of the signal MC sample. At this level of selection, 
we find that the MC models the data well in the elec- 
tron channel of the semileptonic B decay, but less so in 
the muon channel. The disagreement in the muon chan- 
nel appears to derive largely from the continuum back- 
ground and therefore should not affect the real semilep- 
tonic tags. The tag reconstruction efficiency is corrected 
for any data/MC disagreement using a control sample 
described in Section HID. 



B. Selection of B"*" — » t'^u signal candidates 

After the tag B reconstruction, the recoil is studied 

for consistency with the signal modes. All selection cri- 
teria are optimized for each of the different signal r decay 
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modes. The optimization is performed by maximizing the 
signal significance, s/^/s + h, for eacli cliamicl using the 
signal (s) and background (&) MC and assmiiing a total 
branching fraction for i3+ t'^v of 1.0 x lO"**, using 
the PRIM algorithm [14]. This algorithm simultaneously 
optimizes selection criteria over a number of variables by 
relaxing and tightening the constraints on all variables 
until a maximal significance is achieved, allowing only 
up to a fixed percentage of signal and background to be 
removed or restored with each iteration of the selection 
criteria. 

All signal modes contain one charged particle that 
is identified as either an electron, muon, or pion us- 
ing standard particle identification techniques. Both the 
TT^V and the t+ h^ti^V modes contain a pion 
signal track and arc characterized by the number of sig- 
nal tt" mesons. The signal track is required to have at 
least 12 hits in the DCH; its momentum transverse to the 
beam axis, px, is required to be greater than 0.1 GeV/c, 
and its point of closest approach to the interaction point 
must be less than 10.0 cm along the beam axis and less 
than 1.5 cm transverse to it. We demand the invariant 
mass of the signal tt" be between 0.115 and 0.150 GcV/c^. 
The daughter photon candidates must have a minimum 
energy of 50MeV, and their shower shapes are required 
to be consistent with electromagnetic showers. 

Background consists primarily of B^B^ events in 
which the tag B meson has been correctly reconstructed 
and the recoil contains one track and additional particles 
which are not reconstructed by the tracking detectors or 
calorimeter. These events typically contain one or more 
mesons, neutrinos and particles that pass outside of 
the detector acceptance. B'^B^ and continuum events 
contribute background to hadronic r decay modes. In 
addition, some excess events in data, most likely from 
higher-order QED processes (such as two-photon fusion) 
that are not modeled in our MC simulation, are observed. 

Backgrounds are suppressed relative to signal by im- 
posing constraints on the kinematic and shape properties 
of the events. The missing mass is calculated as: 



Mmiss = y (-Br(4S) - -E'vis)^ - (Pr(4S) - Pvis)^. (3) 

Here (i?r(4S)) Pr{4S)) is the four-momentum of the 
T(45), known from the beam energies. The quantities 
E^is and pVis arc the total visible energy and momen- 
tum of the event, which are calculated by adding the en- 
ergy and momenta, respectively, of all the reconstructed 
tracks and photons in the event. Continuum background 
is suppressed with two variables: the cosine of the an- 
gle between the signal candidate and the tag candidate 
thrust vectors (in the CM frame), cos^^, and the mini- 
mum invariant mass constructible from any three tracks 
in an event, M™'". For the background, the cosine of the 
thrust angle peaks at ±1, while the minimum invariant 
mass peaks strongly below 1.5GeV/c^, as can be seen in 
Fig. 2, where the signal and t+t~ background MC are 
shown. We project this 2-d plane into a single variable 



for use in the selection optimization algorithm. The pro- 
jection uses the following empirically derived equation: 

Rcont = ^(3.7 - I cos%|)2 + (M3™/(GeV/c2) - 0.75)2. 

(4) 

Applying selection criteria to i?cont primarily removes 
background from e~^e~ t'^t~, but also suppresses 
other continuum backgrounds. Since the r+ — * 7r+7r°F 
decay proceeds via an intermediate resonance (/9+ 
7r"'"7r'^), further background rejection can be achieved by 
applying requirements on the intermediate meson can- 
didate. In events with more than one recoil 7r°, the 
candidate with invariant mass closest to the nominal tt*^ 
mass [7] is chosen. The invariant mass of the recon- 
structed tt+tt" signal candidates are required to lie be- 
tween 0.64 and 0.86GeV/c2. A quantity analogous to 
cos 6 B-D«£, as defined in section III A, can be calculated 
for r"*" 'n^TT^v by replacing the B with a t and the D^t 
with TT+TT^ in Eq. 2. The analogous quantities of \pt\ and 
Et are calculated assuming the r is from the — > r+i/ 
decay and that the i?"*" is almost at rest in the CM frame. 
We require cos9t—k+-k° < 0.87. 




FIG. 2: (a) Minimum invariant mass of any three tracks and 
(b) |cos^j;| for B+ r+i/ signal MC (solid histogram) and 
e+e^ T^T^ MC (gray histogram). All distributions are 

normalized to niiit area. 



We demand that there are no candidates recon- 
structed in the IFR. For the r+ tt+F channel, we de- 
mand that there arc fewer than two candidate clusters 
in the EMC consistent with being deposited by a K°. 
In the leptonic final states we demand that there are 
two or fewer tt*' candidates. For the t+ e^vV mode, 
we reject events where a photon conversion creates the 
electron by requiring that the invariant mass of the sig- 
nal and tag B lepton pair be greater than O.lGeV/c^. 
We impose mode-dependent selection criteria on the to- 
tal momentum (Pgignai) of the visible decay products of 
the r candidate. 

We further separate signal and background by exploit- 
ing the remaining energy (i^oxtra), calculated by summing 
the CM energy of the neutral clusters (with a minimum 
of 20 MeV in the laboratory frame) and tracks that are 
not associated with cither the tag B or the signal. Signal 
events tend to peak at low i^extra values whereas back- 
ground events, which tend to contain additional sources 
of neutral clusters, are distributed toward higher £^cxtra 
values. The selection applied to E'extra is optimized for 
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the best signal significance, assuming the branching frac- 
tion is 1 X 10^"' and was bhnded for -Eoxtra < 0.5 GeV in 
on-resonancc data until the selection was finalized. 

The signal selection efRciencies for the r decay modes 
are determined from signal MC simulation and summa- 
rized in Table I. The signal efficiencies correspond to 
the fraction of events selected in a specific signal decay 
mode, given that a tag B has been reconstructed. Signal 
selection efficiencies are further corrected by applying the 
factors provided in Table IV which are explained in later 
sections. 

TABLE I: Selection criteria optimized for each signal t decay 
mode. Additional selection criteria are described in the text. 

The signal efficiencies, multiplied by branching fraction, are 
given for each t decay mode, relative to the number of tags. 
Values given in the squared brackets represent lower and up- 
per selection criteria imposed on the respective quantity. 



mode 


e+ 




7r+ 




M„,iss(GeV/c") 


[4.6, 6.7] 


[3.2, 6.1] 


> 1.6 


< 4.6 


Psignal(GeV/c) 


< 1.5 




> 1.6 


> 1-7 


-Rcont 


[2.78, 4.0] 


> 2.74 


> 2.84 


> 2.94 


Sextra(GeV) 


< 0.31 


< 0.26 


< 0.48 


< 0.25 


Efficiency (%) 


4.2 ±0.1 


2.4 ±0.1 


4.9 ±0.1 


1.2 ±0.1 



C. Background Estimation from Bextra Sidebands 

We estimate our backgroimd from the data by study- 
ing events in a sideband region of i?cxtra. We define the 
sideband (sb) region as ii'cxtra > 0.5 GeV, and also define 
signal regions (sig) in E'extra using the appropriate sig- 
nal mode-dependent selection. After applying all other 
selection criteria, we compute from the background MC 
simulation the ratio of events in the sideband (A^MC.sb) 
and signal (A^MC.sig) regions, 



-^MC.sig 
.^MC.sb 



(5) 



Using the number of data events in the sideband 
(-^data,sb) and the ratio R^'~^, the number of expected 
background events in the signal region in data {-^e:xp,sig) 
is estimated, 



expjSig 



ata,sb 



■R 



MC 



(6) 



The sideband background projection (Table II) is taken 
as the number of expected background events. 

The background estimate is validated by performing 
a similar test using sidebands in the mass distribu- 
tion. We select events using D° mass sidebands between 
4(7 and 9<t above and below the nominal D'^ mass, with 
all other signal selection criteria applied. Candidates in 
these regions of the mass distribution are random 
combinations of kaons and pions, and represent a pure 
combinatoric background. We average the yields from 
the upper and lower sidebands and scale this using the 



TABLE II: Comparison of the expected total background, 

computed from data and MC in the D^' mass sideband and 
signal regions, to that computed by projecting the i?extra side- 
band into the Sextra signal region. 





Background Prediction 


signal mode 


e+ n+ 


n+ 




Sextra sideband 


44.3±5.2 39.8±4.4 


120.3±10.2 


17.3±3.3 


D° sideband 


44.2±6.4 42.8±6.0 


113.4±11.6 


16.3±4.5 



ratio of the D° mass sideband and signal region. This 
yields a D'^ mass combinatoric background estimate in 
the D° mass signal region for both data (./V^omb) 
MC (A^^^b)- The remaining component, in the MC, of 
the background which contains real mesons in the tag 
is then computed, 



mMC _ ATlVlU _ AT. 
JV„„„1, — iVt„4.„i — iV, 



peak 



MC 



total 



rMC 



comb ' 



(7) 



and added to the combinatoric component (determined 
from data) to obtain an effective estimate of the total 
background, 



AT-predicted 
-''^ total 



AfMC , j^data 
peak ' comb * 



(8) 



This is done for each reconstructed signal decay chan- 
nel. The method assumes that the background in the 
-Eextra Signal rcgiou can be modeled by the combinatoric 
component of the mass distribution, taken from data, 
and the peaking component of the D° mass distribution, 
taken from MC simulations. Since it uses the mass 
sidebands, it is also statistically independent from the 
-Eextra sideband calculation. 

We find very good agreement between the background 
prediction using the D° mass sidebands and that ob- 
tained from the projection of the -Ecxtra sideband. This 
agreement is demonstrated in Table II and further vali- 
dates our background estimation method. 



D. Correction of tag B yield and iSextra simulation 

The tag B yield and -Boxtra distribution in signal and 
background MC simulation are validated using control 
samples. These samples are further used to define cor- 
rections to efficiencies of selection criteria. "Double- 
tagged" events, for which both of the B mesons are 
reconstructed in tagging modes, B~ D^£~i?gX vs. 
5+ D^t'^v(,X are used as the primary control sam- 
ple. "Single-tagged" events are also used where one B 
decays via B~ D^t~viX and the other B decay is not 
constrained. The double-tagged sample is almost entirely 
free of continuum events. 

We select double-tagged events by requiring that the 
two semileptonic B candidates have opposite charge and 
do not share any particles. We also require that there 
are no additional tracks in the event. If there are more 
than two such independent tag B candidates in the event 
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then the two best candidates are selected as those with 
the largest probabilities of each converging at a common 
origin. The meson invariant mass is shown in Fig. 3 
for the second tag in all double-tagged events. 

We initially determine the tag efficiency using a signal 
MC where one of the two B mesons always decays into 
a generic final state and the other always decays into a 
r+z/ final state. We estimate the correction to the MC 
semileptonic tag efficiency by comparing the number of 
single- and double-tagged events in data and MC. We cal- 
culate the ratio of double-tagged to single-tagged events, 
and we use the ratio of this quantity from data and MC 
as a correction factor for the tag B yield. 

We determine the number of single-tagged events by 
subtracting the combinatoric component under the 
mass peak in events where one B is tagged and the second 
is allowed to decay without constraint (Fig. f). We de- 
termine this component by using Z)° mass sidebands be- 
tween A<j and 7a above and below the nominal D'^ mass. 
A narrower sideband region is used for this study than 
in the background estimate validation due to the com- 
parative flatness of the sidebands in this region and the 
large statistics available at this early stage of selection. 
We then average the yields from these combinatoric 
mass regions and scale by the ratio of the sideband and 
signal region widths. We perform this subtraction using 
events where the meson from one of the semileptonic 
tags is reconstructed as D° K~Tr~^ and the second tag 
decays into any of our allowed final states. The resulting 
single-tagged event yields, and the double-tagged event 
yields, are shown in Table 111. We compared these re- 
sults to that obtained from events where the of at 
least one of the tags decays as D° — > /i'~7r+7r~7r+ and 
found a similar correction. 

We take the uncertainty on the data/MC single-to- 
double-tag ratios as the systematic uncertainty on the tag 
B yield. We find a correction of 1.05 with a 3.6% uncer- 
tainty. This comparison between data and MC provides 
a more realistic environment than signal MC in which to 
compare the various forms of background in the analysis, 
and correct for them. The double-tagged sample alone 
would only correct for B~^B~ backgrounds. 

TABLE III: Single-tag and double-tag yields in data and MC, 
for events where the D° meson from the first tag is required 
to decay as D° — » K~7t~^ . We calculate the ratio of these two 
yields, and take the ratio of these ratios as a correction to the 
tagging efficiency determined from MC. The uncertainty on 
the correction is taken as a systematic error. 



Single Tags Double Tags Ratio 

Data 335417 ±747 1067 ± 33 (3.18 ± 0.10) x 10"^ 
MC 349972 ±572 1065 ± 20 (3.04 ± 0.06) xlO'^ 
Data/MC 1.049 ± 0.038 



We can further test the modeling of E^xtia by com- 
paring it in double-tagged events from data and MC. 




FIG. 3: D° invariant mass from the recoil B meson in double- 
tagged events. On-resonance data (black circles) axe overlaid 
on the combined BB (solid histogram) and continuum (gray 
histogram) MC samples normalized to the data luminosity. 



The -Eextra for the double-tagged sample (Fig. 4) is cal- 
culated by summing the energy of the photons which are 
not associated with either of the tag B candidates. The 
sources of photons contributing to the i^extra distribution 
in double-tagged events are similar to those contributing 
to the -©extra distribution in the signal MC simulation. 

We additionally check the modeling of ii^extra by com- 
paring samples of events where the signal and tag B can- 
didates are of the same sign. We find that for all signal 
modes, there is good agreement between the shape of 
the -Eoxtra from the background prediction and the data 
in this wrong-charge sample. In the pion channel in par- 
ticular, we find the data yield is higher than predicted 
from MC. This suggests that for a pure background sam- 
ple, with a topology similar to that of signal, the -Eextra 
distribution is well-modeled but the background estimate 
cannot be taken directly from the MC background sim- 
ulation. This further validates our choice to take the 
background estimate from the i^extra sideband in data 
and the signal-to-sideband ratio in MC simulation. 



IV. STUDIES OF SYSTEMATIC 
UNCERTAINTIES 

The main sources of uncertainty in the determination 
of the _B+ T^v branching fraction are the tag re- 
construction efficiency (etag)j the efficiency of each sig- 
nal mode (£i), and the number of expected background 
events in the signal region for each signal mode. 

An uncertainty of 1.1% enters the branching fraction 
calculation from the estimation of the number of B^B~ 
events present in the data sample [15]. The tagging ef- 
ficiency and yield in signal MC is corrected using the 
double-tagged and single-tagged samples. The tag B 
yield systematic uncertainty is 3.6%, with a correction 
factor to the yield of 1.05. The systematic uncertainties 
on the signal efficiency depend on the r decay mode and 
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(GeV) 



FIG. 4: -Eoxtra after the reconstruction of two non-overlapping 
semileptonic B candidates. On-resonance data (black circles) 
are overlaid on the combined BB (solid histogram) and con- 
tinuum (gray histogram) MC samples normalized to the data 
luminosity. B"*" — » T+f signal MC (dashed-dotted line) is 
shown for comparison, with arbitrary normalization. 



TABLE IV: Contributions to the systematic uncertainty (in 
percent) on the signal selection efficiencies for different selec- 
tion modes. The total summed uncertainty is added linearly 
with the systematic uncertainties from IFR K" reconstruction 
and -Bextra modeling. The result of this ( "signal B" ) is added 
together in quadrature with the uncertainty on tag B recon- 
struction and the number of BB pairs in the sample (N^^) . 
The "Correction Factor" is a multiplicative factor applied to 
the efficiency for each mode. 



didates is extracted using the double-tagged events, sim- 
ilar to the method used for the i?extra systematic evalua- 
tion. We quantify the data/MC comparison by extract- 
ing the yield with a cut demanding exactly zero (less than 
two) reconstructed IFR (EMC) measured K'l candidates 
remaining, and extracting the yield with a sample where 
any number of K° candidates remain, and take the ratio 
of ratios from the MC and data. The systematic un- 
certainty for vetoing IFR (EMC) candidates is 3.3% 
(3.8%), with a correction factor on the efficiency of 0.99 
(0.97). 

A breakdown of the contributions to the systematic un- 
certainty for each signal mode is given in Table IV. We 
find that the most significant individual effects on the sig- 
nal efficiency arc from the modeling of the -Eoxtra said the 
vetos. The uncertainties on each mode are combined 
by weighting them by the corrected efficiency for a given 
mode, using the efficiencies from Table I multiplied by the 
correction factors given in Table IV. The signal-mode- 
specific systematic uncertainties are summed in quadra- 
ture and then the sum is added linearly with the IFR K'l 
and -Bextra uncertainties, which are correlated among the 
modes. The resulting overall systematic uncertainty on 
the signal efficiency is then added in quadrature with the 
uncertainties on the tag B reconstruction and the num- 
ber of BB pairs in the sample to give a total uncertainty 
of 6.6%. 



V. RESULTS 



T decay mode 
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Total 




6.6 




Correction Factor 


0.951 


0.868 0.964 


0.939 



include effects such as the tracking of charged particles, 
particle identification, and the modeling of tt'^ mesons. 

The systematic uncertainty on the signal efficiency due 
to the mis-modeling of the -Eextra variable is extracted 
using the double-tagged events. We extract the yield 
of candidates satisfying E^extra < 0.5 GeV. This yield is 
then compared to the number of candidates in the full 
sample. Comparing the ratio extracted from MC to that 
extracted from data yields a correction factor, the error 
of which is taken as the systematic uncertainty for Eextra- 
The systematic uncertainty for .Bextra is 3.4% with a cor- 
rection of 0.99. 

The systematic uncertainty on the modeling of can- 



After finalizing the signal selection criteria, we measure 
the yield of events in each channel in the signal region of 
the on-resonance data. Table V lists the number of ob- 
served events in on-resonance data in the signal region, 
together with the expected number of background events 
in the signal region (taken from the Ecxtra sideband pre- 
diction from Table II). Figure 5 shows the Ecxtra distri- 
bution for all data and MC in the signal region, with sig- 
nal MC shown for comparison. Fig ure 6 shows the Ecxtra 
distribution separately for each of the signal modes. 



TABLE V: Observed number of on-resonance data events in 
the signal region are shown, together with number of expected 
background events. 



T 


Expected background 


Observed events 


decay mode 


events 


in on-resonance data 




44.3 ± 5.2 


59 


t"*" — > jJb^VU 


39.8 ± 4.4 


43 


r"*" — > 'k'^V 


120.3 ± 10.2 


125 


T+ ^ 7r+7r°F 


17.3 ± 3.3 


18 


All modes 


221.7 ± 12.7 


245 



We determine the B"*" t^v branching fraction from 
the number of signal candidates si in data for each r 
decay mode, according to Si = NggB{B^ T+z/)£tag£i, 
where is the total number of BB pairs in data. The 
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0.6 0.8 1 

Eext,.(GeV) 

FIG. 5: -Eoxtra distribution. All selection criteria have been 
applied and all signal modes combined. Background MC 
(solid histogram) has been normalized to the luminosity of the 
on-resonance data (black dots), and then additionally scaled 
according to the ratio of predicted background from data and 
MC as presented in section IIIC. -B^ ^ t^v signal MC (dot- 
ted histogram) is normalized to a branching fraction of 10~^ 
and shown for comparison. 




FIG. 



E extra (GeV) 
£extra distribution after all selection criteria for 



_ laj 

e^vu, (b) r+ ^jl^uu, (c) — > it^v, and (d) 
T ' tt'^-k^V. Background MC (solid histogram) has been 
normalized to the luminosity of the on-resonance data (black 
dots), and then additionally scaled according to the ratio of 
predicted background from data and MC as presented in sec- 
tion IIIC. -B"*" T^v signal MC (dotted histogram) is nor- 
malized to a branching fraction of 10~^ and shown for com- 
parison. 



results from each of our four signal decay channels (rich) 
are combined using the estimator Q = C{s + b)/C{b), 
where £{s + b) and C{b) are the likelihood functions for 
signal plus background and background-only hypotheses, 
respectively: 



Lis + fo) ^ n 



- b,r' 



1=1 



(9) 

We include the systematic uncertainties, including those 
of a statistical nature, on the expected background (bi) 
in the likelihood definition by convolving it with a Gaus- 
sian function. The mean of the Gaussian is bi, and the 




Branching Fraction (x 10" ) 




Branching Fraction (x 10 ) 

FIG. 7: (a) Twice the negative natural logarithm of the likeli- 
hood ratio as a function of signal branching fraction hypoth- 
esis and (b) the upper limit as a function of signal branching 
fraction hypothesis (where the horizontal and vertical inter- 
secting lines indicate the 90% CL limit). 



standard deviation (dbj of the Gaussian is the error on 
b^ [16]. 

We calculate the branching fraction central value (in- 
cluding statistical uncertainty and uncertainty from the 
background) by scanning over signal branching frac- 
tion hypotheses between 0.0 and 3.0 x 10^"' in steps of 
0.025 X 10"'' and computing the value of £(s -I- b)/C{b) 
for each hypothesis (Fig. 7a). The branching fraction 
is the hypothesis which minimizes —2 log(£(s -f- b)/£{b)), 
and the statistical uncertainty is determined by finding 
the points on the likelihood scan that occur at one unit 
above the minimum. The systematic error is determined 
as detailed in section IV and computed for the branching 
fraction as a fraction of the central value. 

The upper limit at the 90% CL, including both sta- 
tistical and systematic uncertainties, is determined by 
generating 5000 experiments for each of the aforemen- 
tioned signal branching fraction hypotheses. Each gen- 
erated experiment also includes the expected number of 
background events, and varies the generated number of 
background in each channel according to its uncertainty. 
The total number of events is allowed to vary according 
to Poisson statistics, and systematics are incorporated 
in the efhciency for each channel and the number of B 
mesons originally produced by the collider. The number 
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of signal events in each channel (labeled ?') for each ex- 
periment is thus computed from the branching fraction 
hypothesis as: 

Si=BiX g{ei,Se^) x g(NB+,SNB+), (10) 

where g{x,Sx) represents a number sampled from a 
Gaussian distribution centered on the quantity x with 
systematic uncertainty Sx, and Si and Nb+ are the ef- 
ficiency in each channel and the number of charged B 
mesons produced by the collider, respectively. Each ex- 
periment therefore contains a generated number of signal 
and a generated number of background which will vary 
around the input hypotheses s and b according to the 
above procedures. 

We determine the confidence level of a given signal 
hypothesis by finding the probability that the value of 
the estimator Q in experiments generated according to 
a given composition (signal and background, Qs+b, or 
only background. Qt,) is less than that observed in data 
(Qobs)- The 90% CL limit is determined by using the CLs 
method [17], in which we determine the signal hypothe- 
sis for which P{Qs+h < Qohs)/P{Qh < Qobs) = 1 - 0.90, 
where P(Qs-i-b < Qobs) [PiQb < Qobs)) is the probabil- 
ity that experiments generated assuming a given s + b 
(6) hypothesis have a likelihood ratio lower than that ob- 
served in data (Fig. 7b). 

We determine the branching fraction central value to 

be 

B{B+ T+v) = (0.9 ± 0.6(stat.) ± 0.1(syst.)) x lO""* 

(11) 

and set an upper limit at the 90% CL of 

B{B+ T+v) < 1.7 X 10"^ (12) 

The central value of the branching fraction is in agree- 
ment with that measured by the Belle Collaboration at 
the level of two standard deviations. We interpret this 



result in the context of the Standard Model. Using the 
central value for B{B^ t^v) and taking the known 
values of Gf, tub, "nir and tb [7] we calculate, from 
Eq. 1, /s-IKbl = (7.2l2-O(stat.)±0.2(syst.))xl0-4GeV, 
where the uncertainties are non-Gaussian. Using the 
value of IK&I from [7] we extract /b = 0.167tJ];!]^^ GeV, 
where the uncertainty is dominated by the statistical 
uncertainty on the branching fraction central value. 
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